Maxima of discretely sampled random fields, with an application to ‘bubbles’

A smooth Gaussian random field with zero mean and unit variance is sampled on a discrete lattice, and we are interested in the exceedence probability (P-value) of the maximum in a finite region. If the random field is smooth relative to the mesh size, then the P-value can be well approximated by results for the continuously sampled smooth random field (Adler, 1981; Worsley, 1995a; Taylor and Adler, 2003). If the random field is not smooth, so that adjacent lattice values are nearly independent, then the usual Bonferroni bound is very accurate. The purpose of this paper is to bridge the gap between the two, and derive a simple, accurate upper bound for mesh sizes in between. The result uses a new improved Bonferroni-type bound based on discrete local maxima (DLM). We give an application to the ‘bubbles’ technique for detecting areas of the face used to discriminate fear from happiness. Some key words: Bonferroni; Random fields; Bubbles; Euler characteristic.